Quantum process tomography of a controlled-NOT gate 
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We demonstrate complete characterization of a two-qubit entangling process - a linear optics 
controlled-NOT gate operating with coincident detection - by quantum process tomography. We use 
a maximum-likelihood estimation to convert the experimental data into a physical process matrix. 
The process matrix allows accurate prediction of the operation of the gate for arbitrary input states 
and calculation of gate performance measures such as the average gate fidelity, average purity and 
entangling capability of our gate, which are 0.90, 0.83 and 0.73 respectively. 
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Quantum information science offers the potential for 
major advances such as quantum computing Q and 
quantum communication Q , as well as many other quan- 
tum technologies Q. Two-qubit entangling gates, such as 
the controlled-NOT (cnot), are fundamental elements in 
the archetypal quantum computer 0. A promising pro- 
posal for achieving scalable quantum computing is that 
of Knill, Laflamme and Milburn (KLM), in which linear 
optics and a measurement-induced Kerr-like nonlincar- 
ity can be used to construct cnot gates Gates such 
as these can also be used to prepare the required entan- 
gled resource for optical cluster state quantum computa- 
tion jfj. The nonlinearity upon which the KLM and re- 
lated^, 0] CNOT schemes are built can be used for other 
important quantum information tasks, such as quantum 
nondemolition measurements 0, and preparation of 
novel quantum states (for example, [lflj)- An essential 
step in realizing such advances is the complete charac- 
terization of quantum processes. 

A complete characterization in a particular in- 
put/output state space requires determination of the 
mapping from one to the other. In discrete- variable quan- 
tum information, this map can be represented as a state 
transfer function, expressed in terms of a process matrix 
X- Experimentally, \ is obtained by performing; quan- 
tum process tomography (QPT) [HQ. Q PT has bcen 
performed in a limited number of systems. A onc-qubit 
teleportation circuit an d a controlled-NOT process 
acting on a highly mixed two-qubit state ^3 have been 
investigated in liquid-state NMR. In optical systems, 
where pure qubit states are readily prepared, one-qubit 
processes have been investigated by both ancilla-assisted 
til El and standard QPT. Two-qubit optical QPT 
has been performed on a beamsplitter acting as a Bell- 
state filter [H i. 

We fully characterize a two-qubit entangling gate - 
a cnot gate acting on pure input states - by QPT, 
maximum-likelihood reconstruction, and analysis of the 
resulting process matrix. The maximum likelihood tech- 
nique overcomes the problem that the naive matrix in- 
version procedure in QPT, when performed on real (i.e., 



inherently noisy) experimental data, typically leads to 
an unphysical process matrix. In a previous maximum- 
likelihood QPT experiment E3l ; a reduced set of fitting 
constraints was used. Here we present a fully-constrained 
fitting technique that can be applied to any physical pro- 
cess. After obtaining our physical process matrix, we can 
accurately determine the action of the gate on any arbi- 
trary input state, including the amount of mixture added 
and the change in entanglement. We also evaluate useful 
measures of gate performance. 

The cnot gate we characterize, in which two qubits 
are encoded in the polarization of two single photons, 
is a nondcterministic gate operating with coincident de- 
tectionfthe gate operation presented here is improved 
from PJJ|)- The gate is known to have failed whenever 
one photon is not detected at each of the two gate out- 
puts, and we postselect against these failure modes. This 
gate, described in detail in Ref. ^3; produces output 
states that have high fidelity with the ideal cnot out- 
puts, including highly entan gled states. 

The idea of QPT Q, [H E3 is to determine a com- 
pletely positive map £, which represents the process act- 
ing on an arbitrary input state p: 
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where the A m are a basis for operators acting on p. The 
matrix x completely and uniquely describes the process 
£ , and can be reconstructed from experimental tomo- 
graphic measurements. One performs a set of measure- 
ments (quantum state tomography |2fjj |') on the output of 
an 77-qubit quantum gate, for each of a set of inputs. The 
input states and measurement projectors must each form 
a basis for the set of n-qubit density matrices, requiring 
d 2 = 2 2n elements in each set 0, 0] . For a two-qubit 
gate (d 2 = 16), this requires 256 different settings of in- 
put states and measurement projectors. An alternative 
is ancilla assisted process tomography |l5l flfil l22l| , where 
d 2 separable inputs can be replaced by a suitable single 
input state from a rf 2 -dimcnsional Hilbcrt space. 
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Standard QPT reconstruction techniques typically 
lead to an unphysical process matrix. This is a signif- 
icant problem, as the predictive power of the process 
matrix is questionable if it predicts unphysical gate out- 
put states. For physicality, it is necessary that the map 
be completely positive and not increase the trace. The 
tomographic data can be used to obtain a physical pro- 
cess matrix by finding a positive, Hermitian matrix x 
that is the closest fit in a least-squares sense, and sub- 
ject to a further set of constraints required to make 
sure that \ represents a trace-preserving process (23| : 
E m „ XmnA\A m = I. Practically this is achieved by writ- 
ing a Hermitian parametrization t (Rcf. 20]) of x, and 
minimizing the function 
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where \<f> a ) is the a th input state, \ipb) is the b th mea- 
surement analyzer setting, c a b is the measured number 
of coincident counts for the a th input and b th analyzer 
setting, C is the total number of coincident photon pairs 
within the counting time, A is a weighting factor and 
S is the Kronecker delta. The first sum on the right 
represents a least-squares fit of the Hermitian matrix to 
the data, and the second enforces the set of further con- 
straints when the A k are elements of the Pauli basis. The 
parameter A can be adjusted to ensure that the matrix is 
arbitrarily close to a completely positive map. The tech- 
nique is architecture independent - the photon counts 
can be replaced with the relevant measurement probabil- 
ities for any gate realization. Using our CNOT data, we 
applied a global numerical minimization technique |24j ] 
to find the minimum of /(£) (an alternative maximum- 
likelihood procedure was given in Hj|). 

The ideal CNOT can be written as a coherent sum: 
C/cnot = \{I ® I + I ® X + Z ® I - Z ® X); of ten- 
sor products of Pauli operators {/, X, Y, Z} acting on 
control and target qubits respectively. The Pauli basis 
representation of the ideal CNOT, and our experimental 
process, are shown in Figs.^a) &^b). Physically, the 
process matrix shows the populations of, and coherences 
between, the basis operators making up the gate func- 
tion (note the sign of the coherences corresponds to the 
sign of the terms in C/cnot), analogous to the interpre- 
tation of density matrix elements as populations of, and 
coherences between, basis quantum states. In fact, pro- 
cess matrices are isomorphic with density matrices in a 
higher dimensional Hilbert space [i^, , and the trace- 
preservation condition constrains physical process matri- 
ces to a subspace of physical state density matrices. 



How well does the matrix \ describe the raw data? 
Clearly there will be some discrepancy, as the simple 
matrix inversion (i.e. without maximum-likelihood es- 
timation) produces an unphysical process matrix. It is 
possible to obtain information about the confidence of the 
fit by examining the residuals (Fig. ^c)), i.e. the differ- 
ences A between each of the 256 measurement probabil- 
ities and the corresponding probabilities predicted from 
X- The width of this distribution, er(A) = 0.021, gives an 
idea of the relative error in the process tomography. We 
further test the maximum-likelihood technique by com- 
paring the predicted output state predicted by x) with 
the experimentally determined output state for all 16 in- 
puts. The average fidelity and standard deviation be- 
tween the predicted and measured density matrices are 
0.95 and 0.03, respectively. 

Ultimately, we want to characterize the process relative 
to some ideal: in this case, y™ t, which is the process 
matrix representing C/cnot |28|. We use the process fi- 
delity 0, F> = Tr(xcNOTX), and find F P = 0.87. We 
can obtain a graphical representation of Fp by express- 
ing the process matrix in the "cnot" basis (obtained by 
acting C/cnot on all the Pauli basis elements). In this 
case, Fp is simply the height of the corner (00) element, 
as shown in Fig. ^d). Currently, we arc not able to put 
an error bar on Fp when it is calculated from x- The only 
known technique for obtaining error estimates on the ele- 
ments of the process matrix comes from performing many 
reconstructions of the process matrix with random noise 
added to the raw data in each case ^^|. The incorpora- 
tion of the extra constraints for trace-preservation slows 
down the numerical minimization to the point where it 
is impractical to consider repeating the fitting procedure 
for many simulated data sets. 

The fact that the fidelity of the process is given by the 
height of one element of x in the CNOT basis suggests that 
Fp might be obtained with far fewer experimental set- 
tings than for full QPT. In principle, only d 2 parameters 
are required to find Fp. For our (physically achievable) 
settings j^H , the process fidelity with the ideal CNOT can 
be calculated directly from a 71-element subset of the 
tomographic data. Importantly, any such "direct" rela- 
tionships |29j also allow straightforward error estimates. 
Using this alternative technique, we find F' p = 0.93±0.01 
[30| . The error bar is smaller than F' p — Fp, however, the 
error in Fp is not presently known. 

The average gate fidelity F 0| is defined as the state 
fidelity (3^ between the actual and ideal gate outputs, 
averaged over all input states. There is a simple re- 
lationship between the process fidelity and the average 
fidelity for any process |27j . which we apply with our 
data: T' = (dF' P + l)/(d + 1) = 0.95 ± 0.01. We be- 
lieve that the sub-unit fidelity primarily arises from im- 
perfect mode matching (spatial and spectral overlap of 
the optical beams). Mode mismatch results in imper- 
fect nonclassical interference between control and target 
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FIG. 1: (color) Process matrix of the CNOT gate l3rj|| . (a) Ideal process matrix, in the basis defined by tensor products of 
Pauli operators (imaginary part is identically zero), (b) Maximum-likelihood experimental process matrix, (c) Histogram of 
the differences in probability between experimental data and the maximum-likelihood reconstruction for the 256 tomographic 
measurements. The Gaussian fit 7 exp(— A 2 /a 2 ) has a = 0.021. (d) Real part of the experimentally determined process matrix, 
expressed in the CNOT operator basis, where the 00 element represents the component of C/cnot in the process. The elements' 
imaginary parts are negligible, except for one coherence of magnitude 0.06. In this basis, an ideal CNOT has a unit-valued 00 
element; all other elements are zero. The abscissae are the abscissae of (a) and (b), multiplied by C/cnot- 




photons, and mixture of the individual qubit states as 
well. Mode mismatch is not a fundamental limitation 
for optical quantum gates, and guided mode implemen- 
tations promise an elegant solution. Nonclassical inter- 
ferences with > 99% visibility have been observed with 
single photons and guided mode beamsplitters [3^ . 

Although the fidelity may seem like a simple method 
for comparing processes, it is not ideal, because it does 
not satisfy many of the requirements for a good measure. 
A full list of the desirable properties can be found else- 
where (e.g. Ref. but to some extent they can be 
summarized by the concept that an ideal measure must 
remain valid when used to characterize a gate as part of 
a larger quantum circuit, as well as in isolation. To this 
end, a more appropriate measure, Cp = •v/1 — F' P has 
been developed. Although monotonically related to the 
process fidelity, it has all the properties required. Cp 
is a metric, so that two processes that are identical will 
have Cp = 0; orthogonal processes have Cp = 1. The 
operational interpretation of Cp is that the average prob- 
ability of error Pe for a quantum computer circuit used 
to compute some function obeys Pe < Cp 2 H^j. For our 
gate, P E < 0.07 ±0.01. 

We also introduce a simple relation to characterize how 
much mixture our gate introduces (for details, seeji^l): 
Tr (p 2 ) = (dTr(x 2 )+l)/(d+l) = 0.83, where the quantity 
on the left hand side is the purity of gate output states, 
averaged over all pure inputs. This corresponds to an 
average normalized linear entropy of 0.22. 

An instructive method for examining the action of the 
gate, in terms of fidelity, entanglement (quantified by the 
tangle T 0, |3f|) and entropy, is to make scatter plots 
of these quantities for output states of the gate (Fig. [2J. 
We used ~200000 pure, uniformly distributed (by the 
Haar measure) input states, and the \ matrix, to predict 
a distribution of output states of our experimental gate. 
From this data, we can observe the relationship between 
the mixture added, and both the entanglement gencr- 
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FIG. 2: (a) State fidelity of our CNOT gate outputs (with 
ideal CNOT output states) calculated from \, plotted against 
the linear entropy added by the gate. A perfect experimental 
CNOT process would have F=l, S=0 for all states, (b) Change 
in tangle between input and output, and linear entropy added, 
for our CNOT gate outputs, calculated from \. An ideal CNOT 
would have points with AT distributed between -1 and 1, and 
S=0. 



ated by the gate and the gate output state fidelity. The 
gat e has three separate spatial mode matching conditions 
|l9j, and the contribution of each of these to the overall 
mixture is state dependent, leading to the distribution 
of entropies. There is a clear correlation between the 
fidelity and the amount of mixture added, and the mini- 
mum output state fidelity is 0.83. The shape of the upper 
lobe (AT > 0) of the AT vs. S plot (Fig.^b)) is readily 
understood by the state-dependent mode matching con- 
siderations. States that have the largest change in tangle 
correspond to cases when the gate requires all three mode 
matching conditions to be simultaneously satisfied, and 
since each is not perfectly satisfied, this introduces some 
mixture. There exist partially entangled input states for 
which AT « and all three mode matching conditions 
apply, and these also have higher entropy. When only 
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one mode matching condition applies, the gate cannot 
perform an entangling operation, but only a little mix- 
ture is added. The extension of the lower lobe (AT < 0) 
to AT = — 1 (asymmetric with the upper lobe) can be 
explained by the fact that when the gate acts to disen- 
tangle the input, the addition of mixture also reduces the 
tangle. We find that the maximum increase in tangle of 
the gate (the entangling capability 0|) i s AT max = 0.73. 

In summary, we have demonstrated the full charac- 
terization of a two-qubit entangling quantum process - 
a controlled-NOT gate - by applying physical quantum 
process tomography. With the process matrix, we can 
predict, with approximately 95% fidelity, the action of 
the gate on an arbitrary two-qubit input state. We de- 
termine: an average gate fidelity of 0.90 using the process 
matrix, and 0.95 ± 0.01 using a set of 71 input and mea- 
surement settings; an average error probability bounded 
above by 0.07 ± 0.01; and a maximum increase in tangle 
of 0.73. The main failure mechanism of the gate can be 
observed from the process matrix in the Pauli basis, and 
the scatter plots - some of the operator population is 
incoherently redistributed so that the gate performs the 
identity operation with higher probability than for the 
ideal CNOT, a mechanism that we assign primarily to the 
imperfect mode matching of the interferometers. 
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